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Abstract 



The stability of cosmological solutions in the recently suggested specific mechanism 
of dynamical compensation of vacuum energy is studied. It is found that the solutions in 
^ ' the original version lead to cosmological singularity which could be reached in final (and 

short) time. A modification of the interaction of the compensating field with gravity is 
(N . suggested which allows to escape such singularity. It is shown that generic cosmological 

solution in this model tends to the Friedmann expansion regime even starting from 
initially large vacuum energy. 
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' In a recent paper fl] we considered a mechanism of gross reduction of vacuum energy, 

Q^. Pvac, by a scalar field, <j), coupled to gravity in non- minimal and rather unusual way (see 

' I below Eq. ©). The energy of the condensate of this field in de Sitter background could 

■ diminish p^ac down to the cosmological critical energy density, Pc{t) ~ m?pjt^, transforming 
to '■ exponential cosmological expansion into the usual Friedmann one. The resulting equation 

' of state becomes different from the vacuum one, p = —p. The idea of dynamical adjustment 

of vacuum energy is quite old [5] but no convincing and realistic mechanism had been found 
^ \ for a long time. Though lacking a realistic model, the mechanism of adjustment has two im- 

■ pressive and attractive features: first, the compensation of vacuum energy with the desired 
accuracy of about 100 orders of magnitude and, second, not complete compensation but 
only down to time dependent remnant, ~ m?pjt^. The latter was in fact a generic prediction 
of adjustment models [21 El long before cosmological dark energy had been discovered 1^. 
For a review of the problem of vacuum energy, see e.g. Refs. _5_. 

However, the concrete models of adjustment discussed in the literature [21 El El El suffer 
from numerous shortcomings and realistic cosmology which includes dynamical adjustment 
of vacuum energy has yet to be found. It seems clear that vacuum and dark energies 
are surely related and without solution of the problem of compensation of vacuum energy 
any model of dark energy cannot be considered as complete, though phenomenological 
suggestions [HI can be quite useful. 

In our paper ^ we considered a modification of the model of Ref. 7 and as a result we 
were able to obtain realistic cosmological solutions with dynamical compensation of vacuum 
energy and non-compensated remnants or, better to say, an excessive energy density of (j) 
being cosmological dark energy. To our mind the mechanism described in ^ can be consid- 
ered as moderately successful in the sense that such solutions indeed exist. However, their 
stability has not been studied and its suspected absence may present a serious challenge 
to a possibility of simultaneous solution of both vacuum and dark energy problems in the 
considered frameworks. In this paper we study the stability of the solutions presented in 
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Ref. [T] and find that they are unstable with respect to small perturbations. However, with 
some modification of the underlying Lagrangian governing interaction of (j) with gravity sta- 
bility can be achieved but at the expense of a realistic transition from radiation domination 
stage to matter one. 

The action of the considered model has the form: 



(R + 2A) D,4.D 
2 2i?2 



(1) 



where the metric has the signature (+,—,—,—), is the covariant derivative in this 
metric, and g = — dei[g^y\. We took the units such that TUpi/Sn = 1. Correspondingly the 
cosmological constant A is expressed through the vacuum energy as p^ac = ^mpi/Si: = A. 
The explicit form of the potential C/((/>) is not essential because solutions tend to = 

= const and only the magnitude of the derivative of the potential at this point, U'{(po), 
determines asymptotic behavior of the solution. Hereafter, we take c^o = for simplicity. 

In the cosmological Priedmann-Robertson- Walker (FRW) background equation of mo- 
tion for spatially homogeneous field (j) = takes the form: 

+ u'{4>) = 0. (2) 
while the Einstein equations acquire additional terms related to (p: 




-5m. Pi^) + Pvac] + 2 {g^,,D'- - D^D,) ^^^^ = (3) 

where T^y is the energy-momentum tensor of the usual matter and (Dcp)'^ = Da4>D"(f). 
Taking trace over {p — v) we obtain equation governing evolution of the curvature scalar 

R: 

R-3^-^ + 4 [U{^) + - QD'^-^ = T (4) 

where T = T^. For the spatially homogeneous case the covariant D'Alambertian is = 
/dt'^ + 3Hd/dt and the Hubble parameter, H, is expressed through R as 

R = -6 {2H'^ + if) (5) 

As is found in Ref. these equations allow power law solutions, H = h/t, R = r/t^, 
and [[/((/>) - 

Pvac\ ~ . An interesting feature of this solution is that the term containing 
in Eq. @ is subdominant and can be neglected. However, for the study of stability one, 
as is well known, must retain higher derivative terms in differential equation of motion. We 
checked numerically that the solution of equations ()2I4I5() quickly becomes singular with H 
and/or R going to infinity in finite time. General solution leads to a change of sign of R 
from initially negative value to a positive one and to subsequent collapse [H — > — oo) as 
shown in Fig. ^for = and Fig. I^lfor pm 7^ 0. On the other hand, when the potential 
U is flatter, i.e. smaller U' , the scalar field first changes sign and then all variable {(piH 
and R) blow up for pm = (Fig. 13)).^ In this way a catastrophic cosmological singularity is 
quickly approached. 
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Figure 1: Evolution of the scalar (j), Hubble H and the curvature R for the potential U' = 2.5 
and pm = 0. For initial conditions we take (/)(0) = 0.01, H{0) = 0.01, R{0) = -0.01, 
0'(O) = 10-5 and i?'(0) = 0. 




Figure 2: Same as Fig. ^except for pm(0) = 1. 
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Figure 3: Same as Fig. ^ except for U' = 0.17. 



To cure this unpleasant behavior we changed the Lagrangian in such a way that the 
solution of equations of motion for , -R @ and H © does not lead to the catastrophic 
collapse of the universe for any sign of the curvature R. This may be done if the Lagrangian 
is allowed to depend upon the absolute value of R? If we modify the kinetic term of the 
0-field as 

i?2 ^ R\R\ 

then \H\ would not infinitely rise for any sign of R. Such non-analytical terms in the 
Lagrangian are rather unusual but may be allowed for a toy model. One may have slightly 
better form for non-analytic terms as e.g. y/R"^ + {Dacf))'^- To avoid negative sign of the 
expression under square root the solution should ensure R? > \{D(j))'^\, or if it is not so the 
expression could be properly modified e.g. by changing {Dcf})'^ into (Dcj))^. 
It is convenient to introduce additional notations: 

and w = — — T (7) 



In terms of these new functions and in spatially homogeneous case we obtain the following 
system of equations governing cosmological evolution in the model under consideration: 

■ 0, (8) 



z3 



i + 3Hz - U'{(l)) = 0, (9) 



^ If there exists significant pm the solution behaves as seen in Fig.|5|even for small U' . 
^This modification cannot cure the blow up of (l>,H and R in the case of small U' and pm since this 
instability arises without change of the curvature sign. 
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H + 2H^ + R/6 = 0, 

6D^w + 3^ + M + 4 M [[/(^) + p^^^] + M T = 0. 



(10) 



(11) 



The evolution of the trace of the energy-momentum tensor of the usual matter is determined 
by its equation of state. For relativistic matter T = 0, while for non-relativistic matter: 



The system of equations (|8ll2j) has been solved numerically, see figs. 0] and IHl We 
have found, for more or less general initial conditions, that the curvature, i?, oscillates 
around zero with an increasing frequency and with the amplitude decreasing somewhat 
faster than Possibly particle production by quickly oscillating curvature would make 

the magnitude of R to decay even faster. The Hubble parameter tends to H = l/{2t), as 
is the case of cosmology dominated by relativistic matter. The compensating scalar field 
(p tends to a constant value oscillating around the latter. In the case that non-relativistic 
matter is absent and thus, T = the model corresponds to reasonable Friedmann cosmology 
with the usual expansion law, a{t) ~ t^^^. However, if non-relativistic matter is present and 
T is non-zero it would start to dominate the total matter energy density, T = pm ~ 
while the expansion regime remains relativistic, a ~ 

Still the solutions found in Ref. ^ also exist but to "hit" them a fine-tuning of initial 
conditions is necessary. Moreover, these solutions are unstable and small fluctuations around 
them have one (out of 5) rising mode. It is interesting that if one neglects D^-tevm in Eq. Q), 
as we have done in Ref. P, the "good" solutions found there are stable but the transition 
from relativistic to non-relativistic regime encounters instability. 

So only a moderate success can be reported. The model presented above allows for 
the transformation of de Sitter cosmological solution, driven by vacuum energy, into the 
Friedmann one but the expansion regime is always relativistic with H = l/(2t) indepen- 
dently of the matter content. The matter dominated (MD) solution with H = 2/(3t) is also 
possible but such solution is unstable, small fluctuations around it would drive the system 
back to radiation dominated (RD) regime. Moreover, starting from RD regime, assuming 
that T is small, we have not found any solution which allows for the change to MD regime 
when T begins to dominate matter density. Despite these shortcomings, one may hope that 
the resolution of these problems may be feasible with some modification of the Lagrangian 
describing interaction of the compensating field with gravity. Maybe such constructions 
resembles epicycles in Ptolemaic astronomy but there surely was some truth in the latter. 
Anyhow, adjustment mechanism seems to be the only one that may possibly solve simulta- 
neously both problems of 100 orders of magnitude compensation of vacuum energy and to 
explain the non-compensated remnant in cosmological energy density (dark energy) which 
contributes a relative fraction of order of unity into cosmological energy density at any time 
of the universe history. 
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Figure 4: (a) Evolution of the scalar (f), Hubble h and the curvaturei? for the modified 
kinetic term Eq.®. We take U' = 2.5 and pm = 0. The initial values are i;^(0) = 0.001, 
H{0) = 0.01, R{0) = -0.01, 0'(O) = 10^4 and R'{0) = 0. (b) Evolution of h = Ht and 

U Pvac- 
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Figure 5: Same as Fig. ^except for /Om(0) = 0.1. 
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